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Abstract 

Recent developments of perturbation theory at finite temperature based on effective 
field theory methods are reviewed. These methods allow the contributions from the 
different scales to be separated and the perturbative series to be reorganized. The 
construction of the effective field theory is shown in detail for 4 theory and QCD. It 
is applied to the evaluation of the free energy of QCD at order g 5 and the calculation 
of the g e term is outlined. Implications for the application of perturbative QCD to the 
quark-gluon plasma are also discussed. 
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1 Introduction 



In this paper I review some of the recent developments in perturbative field theory at finite 
temperature that have come about by using effective field theory methods. In particular, 
we will be interested in the behavior of hadronic matter at high temperature. This means 
that the temperature is assumed to be much larger than the mass of the particles involved 
(T ^> m). Another assumption is that the gauge coupling constant g is small (g <C 1); this 
allows us to define three energy scales that satisfy: T ^> gT ^> g 2 T. Hadronic matter is 
expected to undergo a phase transition when T ~ T c ~ 200 MeV between a low temperature 
phase in which it is confined in the form of hadrons and a high temperature phase in which 
quarks and gluons are deconfined. The latter phase, known as Quark-Gluon Plasma (QGP), 
can be studied by perturbative methods since the strong coupling constant a s (T) is expected 
to be small at high temperature. Later on, we will check the validity of this statement; it 
happens that at temperatures of a few times T c , the strong coupling constant is not small 
enough to give rise to a convergent series in the case of the free energy. 

A field theory at high temperature may be described by a lagrangian in which the modes 
with zero Matsubara frequency have been decoupled according to the Appelquist-Carazzone 
theorem Since the resulting theory is 3-dimensional, this method is know as dimensional 
reduction || |[ |J. An efficient approach consists of interpreting the dimensionally reduced 
theory as an effective field theory || whose parameters are to be computed as a perturbative 
expansion in the coupling constant of the original theory. This idea is very attractive when 
studying theories with fermions on the lattice: the dimensionally reduced theory does not 
contain fermions because they have been integrated out completely and this may simplify 
computer simulations |J. The parameters of the effective theory can be determined by 
analytic methods. 

Effective field theory methods have been applied to different problems. Braaten [0 has 
resolved a longstanding problem involving the breakdown of the perturbation expansion for 
the free energy of QCD. Braaten and Nieto have determined the asymptotic behavior of the 
correlator of Polyakov loops || . They have also used this method to carry out explicit calcu- 
lations in 4 || and in QCD [[K], [TTJ] . The problem of the finite T electroweak phase transition 
has essentially been solved in [L2j|: this work has recently been reviewed by Shaposhnikov in 



an article |13[ that is complementary to this review. Andersen [J14] has analytically computed 
the free energy of QED up to order e 5 and the electric screening mass of scalar QED up to 
order e 4 . Arnold and Yaffe [15 have proposed a rigorous nonperturbative definition of the 
Debye screening mass in nonabelian gauge theories. Recently, Karsch et al. [|I^] have com- 
puted the nonperturbative contribution to the g 6 -order term of the free energy of QCD |?|, 
by using lattice simulations. Also, the effective field theory approach has been applied to 
the study of the Minimal Supersymmetric Standard Model . 

The effective field theory approach is now a well-established perturbative method to 
study field theories at finite temperature. It is a useful tool for organizing calculations of 
high-order terms in the perturbation expansion and it is especially powerful in dealing with 
nonabelian gauge theories. It is also an important conceptual tool that explicitly separates 
the different energy scales. In this review I present some of the results of its application to 
QCD. There are two conceptual steps: first, the effective field theory is defined by computing 
its parameters by matching with the original theory; second, the effective theory is used to 
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study problems concerning the high temperature limit of the full theory. 

We will use dimensional reduction to construct an effective field theory of QCD whose 
parameters encode the physics from the scale T; this theory is called Electrostatic QCD 
(EQCD), because it consists of pure Yang-Mills in three dimensions plus a scalar field in the 
adjoint representation which is related to the static mode of the original chromelectric field. 
Since, the effective mass of the scalar is of order gT and magnetostatic fields remain massless, 
we can use the decoupling theorem to construct another effective theory in which the scalar 
field is decoupled. We will end up with a pure Yang-Mills theory in three dimensions that 
is called Magnetostatic QCD (MQCD) because it is only made out of fields related to the 
magnetostatic modes of QCD; its parameters encode the physics from the scale gT. In both 
EQCD and MQCD, the effective lagrangian includes an infinite series of nonrenormalizable 
terms. One of the technically more involved problems in constructing the effective field 
theory is the evaluation of Feynman diagrams. In what follows I will not give details of these 
calculations, I will just state the results of the integrals that are needed and give references, 
which are basically the appendices in [|, [11]]. In all other respects I will try to show the 



construction of the effective theory in detail. 

The next section reviews some concepts that are familiar in the context of quantum 
field theory at finite temperature (more information can be found in the textbooks [jTJ, ID]]), 



introduces the effective field theory approach for 4 theory, and shows a complete calculation 
using its effective field theory. In Sect. |||, the contributions to the QCD free energy from 
the relevant scales are separated. Sect. f| is devoted to the construction of EQCD. The free 
energy of QCD is explicitly computed in Sect. f| up to order g 5 ; also the result at order g 6 
is analyzed. In Sect. [5] we will study the convergence of the perturbative series for the free 
energy. Conclusions and final comments are given in Sect. [7] 

2 Field Theory at Finite Temperature 

The static properties of a system of particles in thermal equilibrium at temperature T are 
obtained from the free energy density 

F = -^logZ, (1) 

where Z is the partition function. All through this paper, we will use the imaginary-time 
formalism which is more suitable than the real-time formalism for studying static properties. 
The partition function of a theory for a field $ is, in general, of the form 



Z{0) = J D$(x,r)exp j-J^dr J dx£($) 



(2) 



where r is the Euclidean time and (3 = 1/T. Also, $ satisfies periodic boundary conditions 
in t if it represents a bosonic field or anti-periodic boundary conditions if it represents a 
fermionic field: 

$(x,0) = ±$(x,/3). (3) 
These conditions allow us to expand $ in its Fourier modes: 

$(x,T)=Tj;^x)e^. (4) 
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The Matsubara frequency of the nth mode is uo n = (2n)7rT for bosons and uo n = (2n + 1)ttT 
for fermions. The term in (|j) with u n = is independent of the Euclidean time r. The 
mode with this zero Matsubara frequency is called static; the rest of the modes involve r 
and are called nonstatic. 

The Feynman rules are the same as in the Euclidean field theory, except that, at finite 
temperature, the time-like component of the momentum is discrete: = (co> n ,k). Conse- 
quently, the loop integration over such a component is replaced by a sum and the rule for 
loops is 



A' 




(5) 



where 3 — 2e is the dimension of the space and u is an arbitrary renormalization scale. 
The factor (e 7 /47r) e is introduced so that, after minimal subtraction of the poles in e due to 
ultraviolet divergences, u coincides with the renormalization scale in the MS renormalization 
scheme. 

For example, in the 4 theory with a massless scalar field and interaction g 2< 3> 4 , described 
by the lagrangian 

c = \m 2 + 9 ^\ (6) 

the leading-order contribution to the self-energy for a particle of momentum P M = (E n , p) 
is given by the tadpole diagram in Fig. |l](a): 



where K 2 = ul + k 2 . 

If we try to compute the next-order correction to the self-energy which is given by the 
Fig. 0(6), we will find out that it is infrared divergent. The same is true for the rest of the 
diagrams shown in Fig. [I]. This mild breakdown of the perturbative expansion may be cured 
by resumming the geometric series of diagrams with one-loop-tadpole insertions like those 
in the Fig. [TJ The result is finite and gives the self-energy at next-to-leading order in g: 

n=«-j , 1 m = mU-^ 9 + o { A. ( 8 ) 

2 Tk K 2 + nW 24 1 4vr vy ; / V ; 



The sum-integral may be found in |20| . We see that the perturbative series is not analytic in 
the coupling constant g 2 , but in g. This is a general feature of the resummation of diagrams 
with insertions of the leading-order contribution of the self-energy. 

Ultraviolet divergences are not a special issue at finite temperature. The very same 
counterterms that would regularize the theory at zero temperature regularize the finite tem- 
perature theory. This can be understood by realizing that a theory at finite temperature 
is defined in a 4-dimensional Euclidean space whose time component is compactified into 
a circle of circumference 1/T. The short- distance behavior is not modified by this sort of 
compactification. 
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2.1 Dimensional Reduction 



Let us consider the Appelquist-Carazzone decoupling theorem for a field theory at zero 
temperature. It states that for a renormalizable field theory with heavy fields (let us say, 
with mass ~ M) and light fields (with mass ~ m <C M), Green's functions with typical 
momentum scale p <C M and that only involve light fields in the external legs may be 
computed without considering heavy field loops. Up to corrections suppressed by powers 
of p/M and m/M, such Green's functions are described by the original lagrangian with the 
heavy fields removed, but with modified values for the coupling constants of the light fields. 
In this sense, the heavy fields decouple from the light fields. 

The partition function of a field theory at finite T may be written in terms of the Fourier 
modes </>„(x), defined in (|), instead of $(x, r): 

Z = J X?0 o (x)X?0„(x)exp j-^dr J dx£(0 o , <t>n, r)| . (9) 

A propagator of the form 1 / '(u^ + k 2 ) can be associated with each mode and the corresponding 
Matsubara frequency represents its mass. In this context, all of the fermionic modes and 
the nonstatic bosonic modes have a mass of order T, while the static bosonic modes are 
massless. 

In the limit of high temperature, nonstatic modes have a large mass of order T when 
compared with static modes, which are massless; therefore, the decoupling theorem suggests 
that they decouple. What remains is an effective theory of the static modes whose partition 
function is of the form 

Z = J X?0(x)exp j-| dx£ cfr (0)j , (10) 



where </>(x) = vT<A>( x ) and £ e fj is the effective lagrangian. The effective theory is defined in 
a 3 dimensional space; this is why this procedure is known as dimensional reduction. Note 
that all fermionic modes have masses of order T and therefore are integrated out. There are 
no effective fields associated with them. Their effects are all incorporated into the effective 
theory through the effective coupling constants for the bosonic modes. 

There is an important difference between the decoupling theorem applied to heavy fields 
at zero temperature and to dimensional reduction at high temperature |3J. For heavy fields 
at zero temperature, the effective theory is usually taken to be renormalizable. It reproduces 
Green's functions of the full theory to all orders in the coupling constant g, up to corrections 
that fall like powers of p/M and m/M. It is possible to reproduce these corrections as well 
but only at the cost of introducing nonrenormalizable terms into the effective lagrangian. 
As we have seen for the <? 2 $ 4 theory at high temperature, perturbative corrections generate 
a mass m for the field of order gT, while the nonstatic field has a mass M of order T. 
The renormalizable effective theory provided by the decoupling theorem reproduces Green's 
functions only up to errors that fall like powers of m/M, but in this case the errors correspond 
to powers of g. Thus the renormalizable theory does not reproduce the full theory to all 
orders in g. Landsman [^T[ has concluded that the validity of the decoupling theorem at 



finite T depends on the theory under study and fails for QCD or 4 theory. This is true 
only if dimensional reduction is taken in the restricted sense of reducing to a renormalizable 
effective theory. As in the case of heavy fields at zero temperature, corrections that fall like 
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powers of m/M can be reproduced by including nonrenormalizable terms in the effective 
theory. 

The effective lagrangian therefore has an infinite number of interaction terms; only a few 
of them are renormalizable, while the rest are non-renormalizable. Each term has a parameter 
that plays the role of an effective coupling constant. These parameters are not arbitrary 
coefficients but instead are completely determined by the condition that the effective theory 
matches the full theory at low momentum. These effective parameters can be calculated in 
powers of the coupling constant g of the full theory. Also, in general, the effective parameters 
depend on an ultraviolet cutoff that cancels the ultraviolet cutoff dependence of the loop 
integrals in the effective theory. In calculations of a physical quantity to a given order in 
g, only a finite number of these parameters will enter. At low orders in g, it is sometimes 
possible to truncate the effective theory to the renormalizable terms or even to the super- 
renormalizable terms. 



2.2 A nontrivial example: the screening mass of 4 theory 

In this subsection we will explicitly construct the effective field theory of a massless field with 
self-interaction g 2 $ 4 described by the lagrangian (|^) and apply it to evaluate the screening 
mass. 

In QED, the electric screening mass describes the asymptotic behavior of the potential 
between two static charges at large distances. The potential is 

f d 3 k e _,kR 

V{R) = J (2^p + n 00 (o,k)' 

where IT^ is the photon self-energy. The asymptotic behavior of the potential is 

-m cl R 

V(R^oo)~— —. (12) 

We observe that the potential at large distances is not a Coulomb potential (1/r) but a 
Yukawa potential (e~ ar /r), so that it goes to zero very rapidly at distances larger that 
l/m e \. That is why m c \ is called screening mass. The asymptotic behavior of the Fourier 
transform in fllTD is determined by the singularity of the integrand closest to the origin, which 
in this case is a pole in the propagator l/{k 2 + IIoo(0,k)); therefore, the electric screening 
mass satisfies |p3| : 

k 2 + n o(0,k) = at k 2 = -m| . (13) 

Analogously, the screening mass m s for the $ 4 theory is given by the location of the pole of 
the static (a? = 0) propagator: 

k 2 + 11(0, k) = at k 2 = -m 2 , (14) 

where U.(u, k) is the self-energy of $. 

The effective lagrangian has to be compatible with the symmetries of the original theory; 
i.e., it has to be symmetric under the exchange <ft ~^ —<P- Then, the most general lagrangian 
is of the form 

£es = \ (V0) 2 + \n 2 </> 2 + -jyA 4 + ^A 6 <f + 5C , (15) 
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where m 2 , A, and \$ are coupling constants of the effective theory. SC represents an infinite 
series of nonrenormalizable terms, each of them with a corresponding effective coupling 
constant. At leading order, any effective parameter will be proportional to T raised to the 
power required by dimensional analysis; the leading power of g 2 is determined by identifying 
the lowest order diagram that contributes to that parameter. The leading-order contributions 
to m 2 , A, and A 6 are shown in Figs. || (a), (b), and (c), respectively. By counting the powers 
of the coupling constant in the diagrams and using dimensional analysis to determine the 
powers of T, their magnitudes must be 



2rp2 



m ~ g 1 
A ~ 
A 6 ~ 



g 2 T 



(16) 
(17) 
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If we are interested in calculations up to order g 4 , we can drop the 6 -term in C e s, since it 
only contributes at higher order. Similarly, the non-renormalizable terms included in 6C can 
also be neglected. What remains is a super-renormalizable effective theory. 

The effective coupling A is computed by matching the full theory and the effective theory. 
Let us consider the action of the full theory defined by the lagrangian (||) 



S 



full 



1/T 



dr / (Fx 



9\i 



2^ 2 + i! $ 



Now, we use the expansion of $ in terms of its Fourier modes 

$(x,r)=T]T0 n (x)e^, 



(19) 



(20) 



and write S'fuU as the sum of two terms; the first of them only depends on the static mode of 
the field and the second term also depends on the nonstatic modes. After integrating over 
r, we get 



full 



drx 



^(V0 O ) 2 + 



+ S S 



full </>0, </>n 



(21) 



Since the first term does not depend on the nonstatic modes, it does not change after 
integrating over the nonstatic modes. As a consequence it can be compared directly with 
the action of the effective theory. After rescaling <fi - 
so that the kinetic terms match, we obtain 



S, 



off 



d x 



T(j) in the effective lagrangian (15) 

(22) 



Comparing ([21]) and (0), we conclude that 

\ = 9 2 T 



(23) 



We have ignored the integration of S\ u n over the nonstatic fields 4> n because we only want to 
obtain A at leading order. This integration contributes to A starting at next-to-leading-order 
and generates the contributions to the other parameters of the effective lagrangian ( |T5] ) such 



as m 2 , A 6 , etc. 
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The actual evaluation of the effective parameters is, in general, performed by matching 
the full and effective theories in the region where they describe the same physics. This region 
corresponds to large distances R ^> 1/T; therefore, we have to match physical quantities 
whose Feynman diagrams have a typical external momentum much smaller than T (p ext <C 
T). In particular, m 2 is determined by matching the electric screening mass m s . In the 
effective theory, the screening mass is defined by 

k 2 + m 2 + n eff (k) = at k 2 = -m 2 , (24) 

where n e g(k) is the self-energy of <fi. We could compute m s in both the full theory and the 
effective theory and then match them to determine m 2 . In the full theory we would have to 
resum infinite series of diagrams such as those described at the beginning of Sect. |2|. 

There is an alternative method that greatly reduces the effort required to determine the 
parameters in the effective theory. It involves introducing an infrared cutoff Am satisfying 
gT <C Am <C T on the internal momenta p- int of diagrams contributing to the physical 
quantity that we want to compute. In the full theory the infrared cutoff should satisfy 
gT <C Am <C T and in the effective theory m <C Am <C A, where A is the ultraviolet 
cutoff of the effective theory. In the full theory with an infrared cutoff p mt > Am, there 
are no complications from infrared divergences since Am ^> gT. We can therefore avoid 
resummations and compute using the strict perturbation expansion in powers of g 2 . In 
the effective theory, with an infrared cutoff p mt > Am, calculations can also be simplified. 
Since the leading order contribution to m 2 is of order g 2 T 2 and Am ^> gT, we can treat 
m 2 as a perturbation and take the propagator of <fi to be simply 1/k 2 . In what follows I 
will refer to the resulting perturbative expansion as the strict perturbation expansion of the 
effective theory. By construction, our effective theory is equivalent to the full theory at low 
momentum. Therefore, the infrared cutoff modifies the full theory and the effective theory 
in precisely the same way. We can therefore match the strict perturbation expansions of the 
full and effective theories to extract the effective parameters. 

Note that we use the strict perturbation expansion only as a device to determine the 
effective parameters. To actually calculate physical quantities we must remove the infrared 
cutoff. In the effective theory, this will require using l/(k 2 + m 2 ) as the propagator of the 
field (f) instead of 1/k 2 . 

Although we have used a momentum cutoff for illustration, we can determine the effective 
parameters by matching the strict perturbation expansions of the full and effective theories 
with any infrared regulator. As we will see, a particularly convenient choice for the infrared 
cutoff is dimensional regularization. While this method does not explicitly cut off the low 
momentum region from internal loops, it still allows the consistent calculation of a strict 
perturbation expansion for both the full theory and the effective theory. 

We now proceed to calculate the mass parameter m 2 in the effective lagrangian ([15]) for 
the $ 4 theory by matching calculations of the screening mass defined by (14) in the full 



theory and by fl24]) in the effective theory. The diagrams in the full theory that contribute 
to 11(0, k) are the ones shown in Fig. |3|. Now, if the g 2 "-order contribution to II(fc 2 ) is called 
H»(fc 2 ), the definition (TTJ) tells us that: 



m 2 



U(k 2 = -m s 2 ) = II« (-m 2 ) + n< 2 > (-m 2 ) + ■ • • . (25) 
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The first term II^ 1 ^— m s 2 ) is given explicitly in (0). It is independent of its argument and of 
order g 2 T 2 . Therefore, at leading order m s ~ g 2 T 2 . Since the argument of lT n )(— m 2 ) is of 
order g 2 T 2 , the strict perturbation expansion for the right side of ( p5| ) is obtained by making 
a Taylor expansion around 0. IT 2 )(— m 2 ) is already of order g 4 , we can set its argument to 
zero up to corrections of order g 4 . We conclude that 



7 2 n 2 



1 

pP 2 



4 Tp P 2 Tp (P 



2\2 



6 r P QP 2 Q 2 (P + Qf 



(26) 



Here and below, the symbol "~" denotes an equality that holds only in the strict perturbation 
expansion. To the order required, renormalization of the coupling constant in the MS scheme 
is accomplished by the substitution 



9_ 

4e V47r 



(27) 



The integrals in (^6|) are evaluated in Ref. || and the final result is 



-g 2 T 2 \l + 

24 y I 



1 , A 

h log 



+2-7+2 



C(-i) 
C(-i) 




(28) 



where A is the mass scale introduced by dimensional regular izat ion. 

The diagrams in the effective theory that contribute to Il e fj(k) are shown in Fig. || (where, 
the propagators and vertices are now those of the effective theory) and in Fig. After 
expanding around k 2 = 0, there is no mass scale in the integrals and the loop diagrams 
vanish in dimensional regularization, 



n eff (0) = Sm 2 . 



(29) 



Consequently, from the definition (J2 



m 2 « 



m 2 + 5m 2 . 



(30) 



We obtain m 2 by matching ( p8|) and (0). Sm 2 in ( |3~0"D is the mass counterterm that 
contains the poles in e that are associated with the mass renormalization. Therefore, it is 
determined to be 

^4 ? ¥W n (31) 
24 y \4irJ e v ; 

The A-dependence from the logarithm in ( PB| ) is partially cancelled by the A-dependence 
of the running coupling constant. The renormalization group equation for the coupling 
constant, 



gives 



9 2 (A) = 9V) 



31 f 

Ait 



) 2log x 



(32) 
(33) 
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which can be used to shift the renormalization scale A to an arbitrary scale fi. We conclude 



™ 2 ( A ) = 7T A 9 2 ^)T 2 1 + 



24 



-3 log -4- + 4 log A 



AtiT 



AixT 



+2-7+2 



C'(-i) 
C(-i) 




(34) 



Note that m 2 (A) depends logarithmically on the ultraviolet cutoff A at order g 4 ; this 
dependence is necessary to cancel logarithmic ultraviolet divergences from loop integrals in 
the effective theory. One should not confuse this ultraviolet cutoff A of the effective theory 
with the infrared cutoff Am. The latter is removed in the matching process. We can explicitly 
see this while using dimensional regularization by introducing two different regularization 
scales A and Am to regularize the ultraviolet and infrared divergences respectively. The 
scale A in fl28|) should be identified with the infrared scale Am. The expression ([30]) for the 
screening mass in the effective theory is then modified to 



m 2 s ^m 2 (A) + -g 2 T 2 



1 (9\ 2 , A 



(35) 



Matching both results the dependence on Am cancels and we recover ( [34]) for m 2 (A). 

Now, having determined m 2 , we can use the effective theory to compute the screening 
mass. Remember that now, the field propagator is l/(k 2 + m 2 ), instead of 1/k 2 . Consider- 
ing (f^j), n e g is given by the diagrams shown in Figs. [3] and |4] evaluated at the point k = im. 
We obtain M 



m 2 {A){ 1 



A 



16nm 



, A 
4 log— + 3 
2m 



8 log 2 



A 



16imi 



(36) 



Now, we can compute the screening mass by inserting (133) and (E3p in 



We find 



+ 



4 log 



4tt\/6 



-l + lllog2-7 + 2 



C'(-i) 
C(-i 




(37) 



where we have set /i = 2ttT. Note that the dependence of m 2 on A in (|3~4[) has been cancelled 
by the ultraviolet cutoff introduced to regularize the effective theory. The term of order g 3 
in (|37|) was first computed by Dolan and Jackiw [22] . The correction of order g 4 was obtained 
by Braaten and Nieto [jj. 

It is worth stressing the distinction between m s and m(A). The screening mass contains 
contributions from both the scales T and gT. In the full theory, it may be obtained by 
resumming the infrared divergent diagrams containing self-energy insertions. On the other 
hand, m(A) only contains contributions from the scale T. By matching (^) and (0) we 
see that it can be obtained by evaluating the screening mass in the full theory by using the 
strict perturbation expansion. We can therefore think of it as the contribution to m s from 
the scale T. 

The technique we have just applied to a massless scalar field with $ 4 self-interaction 
remains essentially unchanged when applied to QCD. In the following sections, we will see 
how it is applied to computing the free energy of QCD. 
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3 Separation of Scales in QCD 



The partition function of QCD is 

Z = J VA„(x,r)VqVqexp j- £ dr j dx£ QCD j , (38) 
where the lagrangian is 

Cqpo^G^G^ + rf^q, (39) 

G a = d^Al — d v A a ^ + gf ahc A h A c v is the field strength, and g is the gauge coupling constant. 
All the quark fields have been assembled into the multi-component spinor q, and the gauge- 
covariant derivative acting on this spinor is = <9 M + igA a ^T a . Quarks are considered to 
be massless; therefore, our description is accurate provided that T is much greater than 
the masses of the quarks considered (m q ). Corrections from the nonzero quark masses are 
suppressed by powers of m q /T. Quarks with masses M q much greater than T can also be 
neglected as they give corrections suppressed by powers of T/M q . 

We will express our calculations in terms of the group-theory factors Ca, Cp, and T F 
defined by 

(40) 
(41) 

(42) 

For an SU(N C ) gauge theory with rif quarks in the fundamental representation with masses 
much smaller that T, these factors are Ca = N c , Cf = (N% — 1)/(2N C ), and Tp = nf/2. The 
dimensions of the adjoint representation and the fermion representation are d,A = — 1 
and d F = N c rif, respectively. 

Now we consider the effective theory that results from integrating out the nonstatic 
modes, which is called Electrostatic QCD (EQCD). Such a theory is only made out of the 
static bosonic modes. The free energy density of QCD, F = —T\ogZ/V can be expressed 
in terms of the effective theory as 



jabc jabd 


= C A 5 cd 


(T a T% 


= CpSij 


tr (T a T b ) 


= T F 5 ab 



„ rpi f log ^EQCD \ 

F = T y* —) ' 



(43) 



where 



^eqcd = J ^ (x)PA(x) exp |- J rfx£ EQCD | . (44) 

and Je is the coefficient of the unit operator that is omitted from the lagrangian. One can 
interpret fpT as the contribution to the free energy from the scale T. 
The lagrangian is 



£eqcd = \G% G% + i(A^)(A^) + \m\A%Al + \\ E {A%A^ + 5£ EQCD , (45) 



where Gr°- = diA^ — djA^+g E f abc A b A'j is the magnetostatic field strength with effective gauge 
coupling constant g E . SCeqcd represents an infinite series of non-renormalizable terms. Note 
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that the effective field theory does not have an effective quark field because fermionic modes 
are integrated out completely; all the effects of the fermions are incorporated into the effective 
parameters. 

Figures |5|(a), (6), (c), and (d) show diagrams that contribute at leading order to the 
parameters of EQCD: fE, m 2 E , g E , and Xe respectively. Counting the powers of g in the 
diagrams and using dimensional analysis to determine the powers of T, we find that the 
magnitudes of the parameters are 

f E ~ T 3 (46) 
m\ ~ g 2 T 2 (47) 
9e ~ 9 2 T (48) 
X E ~ g A T. (49) 

The electrostatic field Aq has a mass of order gT while the magnetostatic fields remain 
massless. Therefore, we can go further in separating the different scales of QCD at high 
temperature by integrating out A . We obtain an effective theory of EQCD which is called 
magnetostatic QCD (MQCD). The free energy density of QCD can be written 

F = T {f E + fM -^f^\ , (50) 



where 



and 



-Zmqcd = J exp j- J dx £ M qcd} • (51) 



mqcd — + &£mqcd • (52) 



1 

4 

The lagrangian is only made out of magnetostatic fields with gauge coupling constant gu\ 
<5£mqcd represents an infinite series of non-renormalizable terms. 

Similarly as we did in the case of EQCD, we can identify the order of the leading con- 
tribution to the MQCD parameters. Figures || (a) and (b) show diagrams in EQCD that 
contribute at leading order to /a/ and g 2 M respectively. Counting the powers of g% in the 
diagrams and using dimensional analysis to determine the powers of m E , we find 

hi ~ m!~s 3 T 3 (53) 
9m ~ 9e~9 2 T. (54) 

We can also see the order at which A# contributes by identifying its leading contribution to 
/at; it is shown in Figure || (c) and turns out to be of order Xettv^ ~ g G T 3 . We see that 
it contributes to the free energy at order g e ; thus, if we are interested in the free energy at 
lower order, we can ignore Xe- Similarly the non-renormalizable terms of EQCD can also be 
omitted. 

The free energy of QCD can be written as 



T 



f E (T,g;A E ) + f M (m 2 E , g E , X E , ■ ■ ■ ] A E , A M ) + f G (g M , ■ ■ ■ ; Am) , (55) 
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where we have defined 

logZ M QCD 

fa = y • (56) 

The factorization scales and Am separate the scales T, gT, and g 2 T. F is written as 
the sum of three terms each of them depending only on the parameters of the corresponding 
theory; their divergences are regulated by the factorization scale parameters. 

We have already seen that the leading contributions to f E and fu are of order T 3 and 
m 3 E ~ g 3 T 3 respectively. The leading contribution to fa can be obtaining by realizing that 
the only parameter with dimensions involved in MQCD is gu and consequently the leading 
contribution is of order g M ~ g 6 T 3 . Since we are interested in computing the free energy of 
QCD up to order g 5 we can ignore the contribution from f G ; then, f E and fu are all that 
we need. We conclude that the free energy of QCD up to order g 5 is given by 



T 



f E {T,g;A E ) + f M (m 2 E ,g E ;A E ) . (57) 



Therefore, we have to determine f E , m 2 E , and g E by matching full QCD and EQCD and then 
calculate f M using EQCD. 



4 The parameters of EQCD 

In this section we will compute the parameters that define EQCD: f E , m 2 E , and g E . They 
are all computed up to the order required to reach our goal of evaluating the free energy of 
QCD at order g 5 . 

4.1 Evaluation of g\ 

In this section, we will compute g E explicitly by comparing the actions of QCD and EQCD. 
To simplify the notation we only consider the Yang-Mills part of the lagrangian (|39|) ; then, 
the action of pure QCD is 

Sqcb = \ fJ T dr J d 3 x G a ^ u G; v , (58) 

where G* = d^A^ — d v A a ^ + gf abc A b ^Al. Now, we use the expansion of A® in terms of its 
Fourier modes, 

A;(x 1 r)=T^(4;)„(x)e- 1 (59) 

n 

and write Sqcd as the sum of two terms; the first of them only depends on the static modes 
of the fields and the second on the nonstatic modes. After integrating over r and rescaling 
(4)o - {A\) /gT, we get 

Sqcd = ^fj d3x G % G % + £qcd[A , A n ] , (60) 

where G?. = di(A°j) - dj(A^) + f abc (A\) Q (Af)Q. Since the first term does not depend on 
the nonstatic modes, it does not change after integrating over the nonstatic modes. As a 
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consequence, it can be compared directly with the action of the effective theory obtained 
from the lagrangian (f45|) . After the rescaling A\ — > (A®) /gE, the action for the effective 
theory is 

Seqcd = J ^ x G%G% ■ (61) 

Comparing (|60[) and (0), we conclude that 

9 2 e = 9 2 T. (62) 

We have ignored the integration of Sqcd over the nonstatic fields {A^) n because we only want 
to obtain g& at leading order. This integration contributes to qe starting at next-to- leading- 
order and generates the contributions to the other parameters of the effective lagrangian ( |45] ) 
such as m%, Xe, etc. 



4.2 Evaluation of 

The effective mass is the contribution to the electric screening mass from the momentum 
scale of order T. In QCD, the electric screening mass m e \ describes the asymptotic behavior 
of the potential between two color charges as in ([12]). The electric screening mass in QCD 
is sensitive to magnetostatic screening effects and therefore requires a nonperturbative 



definition | I5| . However, if one imposes an infrared cutoff that removes those magnetostatic 
effects |§ , one can define a perturbative electric screening mass in terms of the location of the 
pole in the gluon propagator as in ([13]). Although the gluon self-energy is gauge dependent, 



the location of the pole of the propagator is gauge invariant ||24j| . Therefore, if we use a 
gauge-invariant infrared cutoff like dimensional regularization, then the perturbative electric 
screening mass is gauge- invariant. 

In full QCD with an infrared cutoff, the perturbative electric screening mass m e \ is the 
solution to the equation 

k 2 + U(k 2 ) = at k 2 = -m 2 eh (63) 

where U(k 2 ) is obtained from the fi = v = component of the gluon self-energy tensor 
evaluated at k = 0: U^(k = 0,k) = U(k 2 )5 ab . In EQCD with an infrared cutoff, the 
perturbative electric screening mass m e \ gives the location of the pole in the propagator for 
the field Aq(x). Denoting the self-energy function by IlE{k 2 )5 ab , m e \ is the solution to 

k 2 + m| + H E (k 2 ) = at k 2 = -m\. (64) 



By matching the expressions for m e \ obtained by solving ( p3|) and (p4|), we can determine 
the parameter m 2 E . 

We calculate the perturbative electric mass m e \ in the full theory using a strict pertur- 
bation expansion in g 2 and using dimensional regularization with 3 — 2e spatial dimensions 
to cut off both infrared and ultraviolet divergences. Expanding the function iT 1 ^— m 2 ) on 
the right side of (|25f) in powers of m 2 and using the fact that the solution at lowest order 
is m; = nW(0), the resulting expression for the perturbative electric screening mass to 
next-to-leading order in g 2 is 

m 2 el « n«(o) - n« (o)^(o) + n( 2 )(o) . (65) 
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The one-loop diagrams that contribute to U^(k 2 ) are shown in Fig. [7] and two-loop diagrams 
that contribute to H^(k 2 ) are shown in Fig. || their analytical evaluation is detailed in [11 
We find that the strict perturbation expansion for m\ to order g A is 



m 



cl 



g 2 (A)T 2 < Ca + Tp 




. .-> / 5 22 




+ 2 log 



A N 
AtxT 



+ C A T F 3 



4vrT, 

16. n 14 14, A \ 
y log2 + y7 + y log i -^j 



16 



A \ 



+ T^-- T log2-- 7 --log— ) 



6C F T F 




(66) 



In the order g 2 term, we have kept terms of order e for later use. 

The expression (B^) for is an expansion in powers of g 2 . It does not include a g 3 
term, in contrast to the expression for that correctly incorporates the effects of the 



screening of electrostatic gluons [p3 |. This g 3 term arises because the g 4 correction includes 
a linear infrared divergence that is cut off at the scale gT. Since we have used dimensional 
regularization as an infrared cutoff, power infrared divergences such as this linear divergence 
have been set equal to zero. 

In order to match with the expression (|66D, we have to calculate the perturbative electric 
screening mass m e \ in EQCD using the strict expansion in g 2 . Since m 2 E is treated as a 
perturbation parameter of order g 2 , the only scale in the self-energy function IIe(A; 2 ) is k 2 . 
After Taylor expanding in powers of k 2 , there is no scale in the dimensionally regularized 
integrals; therefore, they all vanish. The solution to the equation (|64]) for the perturbative 
electric screening mass is therefore trivial: 



m. 



(67) 



Now, we compare (|66D and (|67|) and take the limit e — >• 0. Note that the expression ( |66|) 
depends on A explicitly through logarithms of A/47rT and implicitly through the coupling 
constant g 2 {A). We shift the scale of the coupling constant from A to an arbitrary scale fi 
by using the renormalization group equation for the running coupling constant 



<? 2 (A) = 9 2 M 



1 + 



2(llCU-47» 



^ ) l0g A 



(68) 



After making this shift in the scale of the coupling constant, the remaining A can be identified 
with the factorization scale A^ that separates the scales T and gT. We conclude that the 
parameter m% is given by 



m| 



g 2 {ti) T 2 \C A + T F 
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^ 2 , 5 22 22 , /x \ 

c ^3 + Y 7 + Y log ;SrJ 



^'.'«|3-^log2 + y 7 + ^lo g ^ 
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+ r^-- T i og 2- 57 --iog^j 



6C F T F 



47T 



(69) 



At this order in g 2 , there is no dependence on the factorization scale A#. 



4.3 Evaluation of Je 



In this subsection, we calculate the coefficient of the unit operator to next-to-next-to- 
leading order in g 2 . The physical interpretation of is that fpT is the contribution to the 
free energy from large momenta of order T. The parameter is determined by calculating 
the free energy as a strict perturbation in g 2 in both full QCD and EQCD, and matching 
the two results. 

In the full theory, the free energy has a diagrammatic expansion that begins with the 
one-loop, two-loop and three-loop diagrams shown in Figs. ||, [K], and [II] Evaluating these 
diagrams in Feynman gauge, we obtain after renormalization of the coupling constant (details 
can be found in |TT| ) 



7l 2 d 




Tj, 



9_ 

47T 



116 220C'(-1 
+ — + 4 7 + 



38 CX-3) ' 

3 C(-l) 3 C(-3), 



146 C(-l) lC'(-3)' 
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c(-i) 3c(-3); 

16 C'(— 1) 8C'(-3)' 



3 C(-l) 3C(-3) 



+ CpTp 



/105 



- 24 log 2^ 



9_ 



(70) 



is a reminder of the strict perturbation expansion of the 



where g = g(A). The symbol ' 
full theory. 

In EQCD, the free energy is given by the expression ([43]). We calculate logZ E q CD using 
the strict perturbation expansion in which g\ and m 2 E are treated as perturbation parameters 
and both infrared and ultraviolet divergences are regularized using dimensional regulariza- 
tion. Since diagrams with massless propagators and with no external legs vanish in dimen- 
sional regularization, the only contribution to logi^QCD which does not vanish comes from 
the counterterm Sfp which cancels ultraviolet divergences proportional to the unit operator. 
The resulting expression for the free energy is simply 



Ue + Sf E ) T . 



(71) 
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The counterterm can be determined by calculating the ultraviolet divergences in logZ E Q CD . 
If we use dimensional regularization together with a minimal subtraction renormalization 
scheme in the effective theory, then SJe is a polynomial in g E , m E , and the other parameters 
in the lagrangian for EQCD. The only combination of parameters that has dimension 3 
and is of order g 4 is g E m E . Thus the leading term in 5Je is proportional to g E m E . The 
coefficient is determined by a 2-loop calculation that is a trivial part of the 3-loop calculation 
in Section 6.1L The result for the counterterm is 



rn AaCa 2 2 1 

6fE = -^f 9EmE ~e 



(72) 



When expressing this counterterm in terms of the parameters g and T of the full theory, 
we must take into account the fact that m E multiplies a pole in e. Thus in addition to 
expression for m E given in ( j69|) , we must also include the terms of order e which can be 
extracted from (|66D. The counterterm (|72|) is therefore 



7T 2 d, 



9 \ rp3 
4:71 , 



Y1C\ - + 2 



«- 1 ' + 21cg. A « 



+ \2C A T F 



^ + l-21og2 + 2^ 



C(-i 
C'(-i 



+ 2 log 




(73) 



Note that minimal subtraction in the effective theory is not equivalent to minimal subtraction 
in the full theory. In addition to the poles in e in (|73|) , there are finite terms that depend on 
the factorization scale A#. 

Matching (|70|) with ( |7T]) and using the expression (fT3~D , we conclude that Je to order g A 

is 




74 (' 

-log2 + 8 7 + - 



3) 



+T F 



20 

— log 

3 4ttT 



C(- 



+ --— log2 + 4 7 
3 5 



■1) 3 C(-3) 
, 16C'(-1) 



8C'(-3) 



3 C(-l) 3C(-3) 



+C F T F 



105 



-24 log 2 




(74) 



where g = g(fi) is the coupling constant in the MS renormalization scheme at the scale fi. 
We have used (|68|) to shift the scale of the running coupling constant from A to an arbitrary 
renormalization scale /i, and we have identified the explicit factors of A that remain with 
the factorization scale A E . 
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5 The Free Energy of QCD 



Once we have understood how to resolve the contributions of the various momentum scales 
in thermal QCD, asymptotic freedom guarantees us that perturbation theory will be under 
control in the high temperature limit. At sufficiently high temperature, the running coupling 
constant will be small enough that calculations to leading order in g will be accurate. How- 
ever, in most practical applications, such as those encountered in heavy ion collisions, the 
temperature is not asymptotically large, and we must worry about higher order corrections. 
The accuracy of the perturbation expansion can only be assessed by carrying out explicit 
perturbative calculations beyond leading order. One of the obstacles to progress in high 
temperature field theory has been that the technology for perturbative calculations was not 
well developed. Only very recently have there been any calculations to a high enough order 
that the running of the coupling constant comes into play. The simplest physical observable 
that can be calculated in perturbation theory is the free energy, which determines all the 
static thermodynamic properties of the system. The running of the coupling constant first 
enters at order g 4 . The free energy for gauge theories at zero temperature but large chemical 
potential was calculated to order g 4 long ago [25]. The first such calculation at high temper- 
ature was the free energy of a scalar field theory with a <fi 4 interaction, which was calculated 
to order g 4 by Frenkel, Saa, and Taylor in 1992 f26fl . (A technical error was later corrected 
by Arnold and Zhai pC§.) The analogous calculations for gauge theories were carried out in 



1994. The free energy for QED was calculated to order e 4 by Coriano and Parwani |27] and 
the free energy for a non-Abelian gauge theory was calculated to order g 4 by Arnold and 
Zhai . The calculation of Arnold and Zhai was completely analytic, and thus represent a 
particularly significant leap in calculational technology. The calculational frontier has since 
been extended to fifth order in the coupling constant by Parwani and Singh and by 
Braaten and Nieto |§ for <p 4 theory, by Parwani [2l| and Andersen for QED, and by 
Kastening and Zhai |3(J and Braaten and Nieto |10, |TT| for non-Abelian gauge theories. In 
the following subsection, the calculation of the free energy for a non-Abelian gauge theory 
to order g 5 based on Ref. [110, O is presented. Also, the calculations that are required to 



obtain the free energy to order g are outlined. 



5.1 QCD free energy: up to order g h 

In order to reach our goal of obtaining the free energy of QCD to order g 5 , it remains to 
evaluate the coefficient of the unit operator of MQCD, Jm, which is the contribution of the 
scale gT to the free energy. 

Through order g 5 , ju is proportional to the logarithm of the partition function for EQCD: 

log^EQCD 

/m = y • (75) 

The lagrangian of EQCD that we are considering at this point is 

£ E qcd = -G%G% + \{D t A,f + \m\Al , (76) 

where now the parameters Qe and rriE are known. In order to calculate the contribution to 
fu using perturbation theory, we must incorporate the terms in the lagrangian that provide 
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electrostatic screening into the free part of the lagrangian. The necessary screening effects 
are provided by the AqAq term in the EQCD lagrangian. Thus we must include the effects 
of the mass parameter m E to all orders, while treating all the other coupling constants of 
EQCD as perturbation parameters. The only coupling constant that is required to obtain 
the free energy to order g 5 is the gauge coupling constant g E . 

The contributions to log ^eqcd of orders g 3 
2- loop, and 3-loop diagrams in Figs. |T2 
the propagators of the A 



g A , and g 5 are given by the sum of the 1-loop, 



[E| and |TJ. The solid 
fields 



wavy, dashed lines represent 
and the associated ghosts, respectively. We 



field, the Ai 

evaluate these diagrams in Feynman gauge. The details of this calculation may be found 
in where the methods developed by Broadhurst were used to evaluate analytically 



the integrals involved by the 3-loop diagrams. The resulting expression for the logarithm of 
the partition function is 



dt 
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(77) 



where 8f E is the counterterm associated with the unit operator of the EQCD lagrangian and 
A is the scale of dimensional regularization. It can be identified with the ultraviolet cutoff 
of EQCD. The ultraviolet pole in e in the term proportional to g E m E in ( [77] ) is cancelled 
by the counterterm 8f E , which is given in (|72|). Our final result is therefore 



/jw(Ai 



3(4tt) 
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(7f 



The coefficient J'm in (f78|) can be expanded in powers of g by setting g E = g 2 T and 
by substituting the expression ( |69"D for m E . The complete free energy to order g 5 is then 
F = (/e + fu)T. Note that the dependence on the arbitrary factorization scale A# cancels 
between f E and Jm, up to corrections that are higher order in g, leaving a logarithm of 
T/m E . This g A \og(g) term is associated with the renormalization of f E , and its coefficient 



can be determined from the evolution equation [11 



A 



d 



E 



dA, 



d A C_ 



A 2 



g E m E + 0{g^ 



(4tt) : 



(79) 



There is no g 5 log(g) term in the perturbation expansion for F, and this is a consequence of 
the vanishing of the g E term in the evolution equation for m\ analogous to (|79|). 



5.2 QCD free energy: beyond order g 5 

In 1980, Linde pointed out a problem concerning perturbation theory at finite temper- 
ature. Let us consider diagrams that contribute to the free energy of QCD like the one in 
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. L+l 2L i 

'-^in^n (so 



Fig. [15] with L + l loops. Its infrared behavior is described by an integral of the form 

L+l 2L 

, , k 2 + m 2 

1=1 j=l 3 

where we have inserted a mass m into the propagator as an infrared cutoff. From dimensional 
analysis we see that when 

• L < 3, I is infrared finite. 

• L = 3, I ~ # 6 T 4 log(T/m). 

• L > 3, / ~ g 6 T A (g 2 T/m) L - 3 . 

This means that if we consider nonstatic gluons whose mass m is set by the Matsubara 
frequency which is of order T, then I ~ g 6 T 4 (g 2 ) L ~ 3 ; the magnitude of the contribution 
decreases by a factor g 2 for each loop as in ordinary perturbation theory. If we consider 
electrostatic gluons whose screening mass is of order gT, then / ~ g 6 T 4 g L - 3 ; the magnitude 
of the contribution decreases by a factor of g for each additional loop. However, if we consider 
magnetostatic gluons to be screened at a scale m of order g 2 T, I ~ g 6 T 4 ; all of the diagrams 
contribute to order g G , no matter how large the number of loops is. 

This problem remained open for many years. In 1994, Braaten |7J proposed a solution of 
this puzzle in the context of the effective field theory approach we are reviewing here. Later, 
Braaten and Nieto [FJ] analyzed the diagrams that contribute to order g 6 from the different 



scales. The free energy is given by (|55|) 

F = T(f E + f M + f G ) , (81) 

where f E , fni, and fa give the contributions to the free energy from the scales T, gT, and 
g 2 T respectively. We proceed to outline the calculations that would be required to obtain 
F to an accuracy of g 6 . In the full theory, we have to calculate contributions to f E from 
4-loop diagrams and, also, the terms up to order g 4 and eg 4 for g 2 E and A^- In EQCD, we 
need to calculate the 4-loop diagrams that give the term g E in f M and the 1-loop diagram 
that gives the \ E m E term. Also, g E has to be computed up to order g 4 . Finally, there is a 
contribution from MQCD which can be written in the form 

fa= (a + b\og^-)g 6 M . (82) 

The number b may be calculated by evaluating the logarithmic ultraviolet divergence in 4- 
loop diagrams of MQCD, but the number a requires nonperturbative calculations. It can 
be calculated using lattice simulations of pure gauge theory in three dimensions. Recently, 



Karsch et al. [[U| have reported lattice calculations that estimate the value of a. 

The contributions to the free energy at higher order in g can be analyzed in a similar 
way. The contributions from the scales T and gT can all be obtained from diagrammatic 
calculations in full QCD and EQCD. The contributions from the scale g 2 T require nonper- 
turbative calculation in MQCD. It is clear that one can write an expansion for F in powers 
of g at arbitrary order. In this sense, there is no breakdown of the weak-coupling expansion 
at any order. 
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6 Convergence of Perturbation Theory 



We have calculated the free energy as a perturbation expansion in powers of g to order g 5 

F=(f E + f M )T, (83) 



where Je is given by (|74| ) and J'm by (|78|) . In this section, we examine the convergence of 
that perturbation expansion; the analysis is based on Refs. [0, [TTH . For simplicity we focus 



on the case of QCD with nf flavors of quarks. The expansion of the free energy in powers 
of ^JoT s with a s = g 2 / is 
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(89) 



The coefficient F 2 was first given by Shuryak |]33| . The coefficient of F3 was calculated 
by Kapusta |34] and Kalashnikov and Klimov |35| . The coefficient of order a 2 log a s was 
obtained by Toimela |36| and the coefficient of order a 2 by Arnold and Zhai [f20" |. The 
coefficient F5 has also been calculated independently by Kastening and Zhai p0|. 



We now ask if the expansion (Q) is well-behaved. If the series is apparently convergent, 
then it can plausibly be used to evaluate the free energy. We study the expression ( |54"[ ) 
at different temperatures T > T c ~ 200 MeV. We choose the renormalization scale to 
be /i = 2ttT, which is the mass of the lightest nonstatic mode. In Table [I] we give the 



contributions coming from each order in 



'a. 



rescaled so that the leading order term is 1. 



We see that the correction of order cfij 2 is the largest unless the temperature T is greater 
than 2 GeV. The term of order a z J 2 is smaller than the term of order a s only when the 
temperature is greater than about 1 TeV. 
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T (GeV) a s (27rT) expansion for F 

O250 (L321 1 - 0.282 + 0.583 + 0.276 - 1.094 

0.500 0.239 1 - 0.210 + 0.374 + 0.010 - 0.524 

1 0.194 1 - 0.167 + 0.267 + 0.033 - 0.287 

2 0.165 1 - 0.142 + 0.209 + 0.011 - 0.191 
1000 0.074 1 - 0.0624 + 0.0619 + 0.0106- 0.0242 

Table 1: Perturbation expansion for the Free Energy F in units of (— 8ir 2 T A / 45) Fq at different 
temperatures. 



T (GeV) a s (27rT) expansion for f E 

0.250 0.321 1 - 0.282 + 0.489 

0.500 0.239 1 - 0.210 + 0.271 

1 0.194 1 - 0.167 + 0.132 

2 0.165 1 -0.142 + 0.096 

Table 2: Perturbation expansion for f E in units of (— 87r 2 T 4 /45)F at different temperatures. 



We can go further in our analysis provided that we have separated the contributions from 
the scales T and gT. We proceed to study the perturbation expansion at the scale T. The 
term f E which gives the contribution to the free energy from the scale T is given in (ff4|): 
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(90) 



Again we choose the renormalization scale to be \x = 2ttT. We chose the factorization scale 
to be 27rT, to avoid large logarithms of A e /2ttT. Rescaling the contributions from each 
order in a s , we obtain the expansion of f E shown in Table 0. We see that the term of order 
a s is reasonably small for all the values of T that are given. Note that one could choose 
so that the contribution of order a 2 s cancels. Therefore, the size of the order a 2 s correction is 
not a very good test of the perturbation expansion. 

It is interesting to study the convergence of the other parameter of EQCD that we have 
computed, m E , which is given by (|6 



4tt a s (fi) T 2 \ 1 + |n/ + 



0.612 - 0.488n/ 



0.0428n^ 
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T (GeV) a s {2itT) expansion for m E 

(L250 (L321 1 - 0.124 

0.500 0.239 1 - 0.093 

1 0.194 1 - 0.098 

2 0.165 1 - 0.083 

Table 3: Perturbation expansion for m 2 E in units of 4iia s (27iT)T 2 (l +nf/6) at different 
temperatures. 



T (GeV) a s (27rT) expansion for / 



M 



0.250 0.321 1 - 1.050 - 1.694 

0.500 0.239 1 - 1.047- 1.262 

1 0.194 1 - 0.947- 0.930 

2 0.165 1 - 0.935 - 0.791 



Table 4: Perturbation expansion for f M in units of -2/(37r)[47ra,(27rT)T 2 (l + n f /6)} 3 / 2 at 
different temperatures. 
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Again, setting \x = 2ttT we obtain the corrections to the leading order result that are sum- 
marized in Table Here, we also see that the next-to-leading order correction is reasonably 
small for all the values of T. Based on these results, we conclude that the perturbation series 
for the parameters of EQCD are well-behaved even at temperatures as low as 250 MeV. 

We next examine the behavior of the perturbation expansion for EQCD. The term Jm is 
given by (|7F 



/iw(Ai 
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We choose the renormalization scale of EQCD to be uie, which is mass of the electrostatic 
mode. In Table § we give the contribution at each order in g E /rriE for the factor in square 
brackets in fl9~2|). The next-to- leading order correction could be made to vanish by a suitable 
choice of A#. Therefore, it is not a good test of the convergence of the expansion. The next- 
to- next-to-leading order correction is independent of A^ and is smaller than the leading order 
term only if the temperature T > 2 GeV. Thus, the temperature for which the perturbation 
series for ju is well-behaved is much higher than that required for the parameters of EQCD 
to have well-behaved perturbation series. 

This analysis suggests that the slow convergence of the expansion for F in powers of y/ctl 
may be attributed to the slow convergence of perturbation theory at the scale gT. 
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7 Conclusions 



In this paper, we have reviewed effective-field-theory methods to study the high T limit of 
QCD. These methods have been used to unravel the contributions to the free energy of QCD 
at high temperature from the scales T, gT, and g 2 T. Also, the free energy has been explicitly 
computed to order g 5 and the calculation of the g 6 contribution outlined. The calculation 
was significantly streamlined by using effective-field-theory methods to reduce every step of 
the calculation to one that involves only a single momentum scale. 

Our explicit calculations allow us to study the convergence of the perturbation expansion 
for thermal QCD. They suggest that perturbation theory at the scale gT requires a much 
smaller value of the coupling constant than perturbation theory at the scale T. At the scale 
T, perturbative corrections are small for all temperatures T > T c ~ 200 MeV. Of course, 
even if this condition is satisfied, the perturbation expansion may break down anyway, but 
this is certainly a necessary condition. At the scale gT, perturbative corrections can be small 
only if T > 2 GeV. Thus, in order to achieve a given relative accuracy, the temperature T 
must be an order of magnitude larger for perturbation theory at the scale gT compared to 
perturbation theory at the scale T. 

There is a range of temperatures in which perturbation theory at the scale gT has broken 
down, but perturbation theory at the scale T is reasonably accurate. In this case, one can still 
use perturbation theory at the scale T to calculate the parameters in the EQCD lagrangian. 
However, nonperturbative methods, such as lattice simulations of EQCD, are required to 
calculate the effects of the smaller momentum scales gT and g 2 T. While one could simply 
treat the entire problem nonperturbatively using lattice simulations of full QCD, the effective- 
field-theory approach provides a dramatic savings in resources for numerical computation. 
The savings come from two sources. One is the reduction of the problem from a 4-dimensional 
field theory to a 3-dimensional field theory. The other source of savings is that quarks are 
integrated out of the theory, which reduces it to a purely bosonic problem. 
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Figure Captions 

1. Diagrams that contribute to the self-energy of 4 . 

2. Leading order contributions to the effective parameters of £ e g for <f) 4 . 

3. Diagrams that contribute to the self-energy up to 2- loop order for 4 . 

4. Diagrams that contribute to the self-energy up to 2-loop order for 4 involving the 
mass counterterm. 

5. Leading order contributions to the parameters of EQCD: (a) f E , (6) m E , (c) g E , and 
(d)X E . 

6. Leading order contributions to the parameters of MQCD: (a) /m, (b) g 2 Ml and (c) 
the first contribution to J'm that involves X E . Solid and wavy lines represent the 
propagators of the A Q field and the Ai fields, respectively. The solid blob represents 
the vertex associated with \ E . 

7. One-loop Feynman diagrams for the gluon self-energy. Curly lines, solid lines, and 
dashed lines represent the propagators of gluons, quarks, and ghosts, respectively. 

8. Two-loop Feynman diagrams for the gluon self-energy. The solid blob represents the 
sum of the one-loop gluon self-energy diagrams shown in Fig. [f[ 

9. One-loop Feynman diagrams for the free energy of QCD. 

10. Two-loop Feynman diagrams for the free energy of QCD. 

11. Three-loop Feynman diagrams for the free energy of QCD. 

12. One-loop Feynman diagram for the logarithm of the partition function of EQCD. 

13. Two-loop Feynman diagram for the logarithm of the partition function of EQCD. 

14. Three-loop Feynman diagrams for the logarithm of the partition function of EQCD. 

15. Diagram that may give rise to the breakdown of perturbation theory in thermal QCD. 



26 



(b) 




Fig. 1 




Fig. 2 




Fig. 4 












Fig . 5 
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